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Abstract. Recent analyses of cosmic ray arrival directions have resulted in evidence for a 
positive correlation with active galactic nuclei positions that has weak significance against an 
isotropic source distribution. In this paper, we explore the sample size needed to measure a 
highly statistically significant correlation to a parent source catalogue. We compare several 
scenarios for the directional scattering of ultra-high energy cosmic rays given our current 
knowledge of the galactic and intergalactic magnetic fields. We find significant correlations 
are possible for a sample of >1000 cosmic ray protons with energies above 60 EeV. 
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1 Introduction 

The origin of ultra-high energy cosmic rays (UHECRs) is a long-standing unsolved problem, 
which has defied an observational solution in large part due to magnetic field scattering. As 
cosmic rays propagate towards Earth, the galactic and intergalactic magnetic fields deflect 
their trajectories so that only the highest energy particles are rigid enough to point back to 
their sources. UHECRs also lose energy as they propagate through space due to interactions 
with photon backgrounds [1, 2], meaning that the sources of UHECRs observed with ener¬ 
gies ^ 6 X 10^® eV are expected to be nearby. The intergalactic magnetic fields, although 
weak, interact with UHECRs over the whole course of their trajectory, which can result in 
large deflections. These combined effects mean that arrival directions can become isotropized 
resulting in weakly significant correlations with a parent source catalogue. 

Recently, the Auger collaboration has been searching for correlation in the arrival di¬ 
rections of ultra-high energy cosmic rays with the active galactic nuclei (AGN) distribution 
represented by the Veron-Cetty & Veron (VCV) catalogue [3-6]. They compare their cat¬ 
alogue correlation with the correlation due to an isotropic source distribution to estimate 
the statistical significance of their result. The test assumes three parameters: 1) the energy 
threshold, motivated by the expectation that above the Greisen-Zatsepin-Kuzmin (GZK) cut¬ 
off [1, 2] UHECRs are extra-galactic, 2) a correlation angle, motivated by the degree to which 
the galactic magnetic field may scatter events and result in significant source confusion, and 


- 1 - 




3) the maximum redshift of catalogue sources to cousider iu the correlatiou aualysis, mo¬ 
tivated by the expectatiou that UHECRs above the GZK cutoff are withiu a few huudred 
megaparsecs. 

Ill their latest release [6], the correlatiou to VCV AGN provided ouly weak evideuce of 
auisotropy aloug with weak clusteriug arouud Geutaurus A. The HiRes aud Telescope Array 
collaboratious performed similar studies fiudiug uo highly statistically siguificaut correlatiou 
with a source catalogue [7-9]. 

Future UHEGR detectors could provide the exposure ueeded to reveal auisotropy aud 
clusteriug to source catalogues. Proposed space-based observatories such as the JEM-EUSO 
missiou [10] aud radio detectiou iustrumeuts [11-13] offer the possibility of siguificautly ex- 
teudiug the sample of UHEGRs available for source correlatiou aualysis. 

This paper aims to estimate the UHEGR sample size that could lead to a statistically 
siguificaut source catalogue correlatiou with a full-sky survey of UHEGRs above the euergy 
threshold of 60 EeV. Giveu current limitatious ou knowledge of source composition as well 
as galactic and intergalactic magnetic fields, we provide estimates for various assumptions on 
these parameters. 

There have been several simulation studies on the effects of cosmic ray scattering by 
galactic and intergalactic magnetic fields and their impact on source identification efforts [14- 
21]. The approach presented here is not to provide another detailed simulation of cosmic ray 
propagation. Instead, we use existing simulation results to provide a parametric simulation 
that readily ties the behavior of composition, galactic and intergalactic magnetic field param¬ 
eters to the catalogue correlation analysis presented by the Auger collaboration. 

A related study explored the sensitivity of a JEM-EUSO-like instrument to detecting 
anisotropic arrival directions of UHEGRs, assuming several different astrophysical source 
distributions, but a single magnetic field scattering model. They report that anisotropy 
would be observed in the lighter component of the UHEGR population with a detector that 
has an order of magnitude increase in exposure over current experiments [22]. Our work 
complements this approach by considering several magnetic field models, but with a single 
presumed source distribution. 

There have been other studies for all-sky instruments that employ different techniques for 
identifying anisotropy. One study employs the two-point correlation function to estimate the 
sensitivity of detecting multiplets, using the 2MRS catalogue [21]. Another study estimates 
the sensitivity to excesses of the dipole and quadrupole moments in the arrival direction 
distributions above an isotropic distribution [23]. These techniques do not employ the use 
of a source correlation catalogue to estimate the degree of anisotropy. In contrast, our work 
focuses on estimating the sensitivity to detecting a statistically significant correlation between 
UHEGR arrival directions and a source catalogue, relative to an isotropic distribution. In the 
best possible case, the source catalogue and the correlation catalogue are identical. 

The paper is organized as follows. In section 2 we cover UHEGR scattering due to 
galactic and extragalactic magnetic fields. Section 3 presents a parametric model for cosmic 
ray propagation. Section 4 presents simulated correlation results of UHEGRs with a source 
catalogue under different assumptions of galactic and extra-galactic magnetic fields. 

2 Cosmic Ray Scattering due to Galactic and Intergalactic Magnetic Fields 

The magnetic fields affecting cosmic ray deflection are not entirely understood. The inter¬ 
galactic magnetic field (IGMF) has a wide range of uncertainty in its parameters given how 
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difficult it is to observe its effects. The parameters that describe the IGMF are the magnetic 
held strength Bq and coherence length \b- A long coherence length \b means the magnetic 
held is constant in direction and magnitude over large distances, while a short \b indicates a 
turbulent magnetic held that is varying direction and magnitude on small scales. Depending 
on the distance of propagation, cosmic ray scattering can depend as much on the magnetic 
held strength as it does on coherence length. The galactic magnetic held (GMF) has turbulent 
contributions from the disk and the halo, which are also poorly constrained. In the follow¬ 
ing we discuss the parameters of the intergalactic and galactic magnetic held and describe a 
parametric model of their scattering of ultra-high energy cosmic rays. 

2.1 Intergalactic magnetic fields 

Gurrent bounds in the strength of the intergalactic magnetic held constrain Bq to a range 
10“^^ —10“® G. The upper bound Bq < 10“® G is due to the impact of intergalactic magnetic 
helds on cosmological perturbations and GMB anisotropies using Planck data [24]. The lower 
bound Bq > 10“^'^ G is due to the non-observation of GeV y-rays by the Fermi Large Area 
Telescope following from TeV y-rays observed by HESS [25]. 

The intergalactic magnetic held coherence length is also poorly constrained over a 
large range. A theoretical argument of magneto-hydrodynamic turbulence decay results in a 
constraint that Xb > 0.1 Mpc at Hq = 1 nG and Xb > 10“® Mpc at B = 10“^® G. A detailed 
review of constraints on the IGMF parameters can be found in [26]. 

A recent study [27] reports evidence of a redshift dependence on the rotation measure 
from quasar light curves. The observations result in a measurement of intergalactic magnetic 
held strength Bo = l]']g 3 nG. However, similar works did not hnd a strong correlation between 
rotation measure and redshift, due to the large variation in the intrinsic rotation measure of 
each source [28-30]. This means that Bq could be well below 1 nG. As will be shown later, 
the IGMF, even at this strength, is not the dominant source of scattering. 

Assuming a Bq = 1 nG combined with the MHD constraints [26], we arrive at a lower 
bound on the coherence length Xb >0.1 Mpc. It is interesting to note that while Xb = 0.1 
Mpc could be generated by AGN winds [26, 31], longer coherence lengths can currently only 
be explained by a cosmological origin of primordial magnetic helds [26]. 

2.2 Galactic magnetic fields 

The galactic magnetic helds, in addition to being described in terms of their strength Bq 
and coherence length Xb, are also bound over a distance D. The galactic magnetic helds are 
classihed in terms of the disk and halo contribution, each with their own parameters Bq, Xb, 
and D. The disk and halo magnetic helds can have regular (D < A^) and turbulent (D > Xb) 
components. 

The galactic disk’s regular magnetic held effect on the scattering of UHEGRs was studied 
by Stanev [32] with a magnetic held strength of Bq ~2 ^uG. The distance over which the 
UHEGR is dehected is limited by the thickness of the disk and assumed to be D ~2 kpc. The 
turbulent component of the galactic disk’s magnetic held is assumed to have a strength of 
Bo ~4 /iG with a coherence length Xb ~50 pc [33]. The magnetic held of the galactic halo is 
less well known. Studies conducted by Jansson et al. 2009, [34] estimate a regular magnetic 
held with strength Bo ~2 /rG over a distance of D ~8.7 kpc. Other measurements [35] indicate 
the halo magnetic held strength may be signihcantly higher. No observational constraints on 
the turbulent component of the galactic halo magnetic held are known. 
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2.3 Scattering of UHECRs by Magnetic Fields 


Several parameterizations of the magnetic field scattering of UHECRs exist in the literature 
[36-38]. In this work, we adopt the results from Lee et ah, 1995 [39], as presented in Neronov 
and Semikoz, 2009 [33], which provide a parameterization valid for energies E > 10 EeV, 
including the varying scales for the coherence length Xb, magnetic held strength Bq, and 
charge number Z. This parameterization for the mean scattering angle of a UHECR due to 
interactions with a magnetic held is given by 
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for a turbulent held. 

UHECRs are scattered over a distance, D, which varies with magnetic held model. For 
intergalactic magnetic helds, D is the distance to the source, while for galactic magnetic helds, 
it is the bounding distance relevant to the different regions of the galaxy. Here we average 
over the structure in the galactic halo and disk. A more detailed treatment is described in 
[39]. 

It is important to note that this parameterization is valid for small dehection angles 
{'&scat ^ 10°). We do not include a separate parameterization for higher scattering angles 
since it is not relevant to this study. As will be shown later in this work (see Sec. 4.2), only 
dehections well within the range of validity of this parameterization are shown to correlate 
to their sources. 

Table 1 summarizes the magnetic held scales considered here and includes a represen¬ 
tative scattering angle for a proton with energy E = 10^° eV. The energies used for this 
parameterization are the energy of the UHEGR upon entering the magnetic held. For IGMF 
propagation, this is the energy of the UHEGR at the source, while for galactic magnetic held 
scattering E is the energy of the UHEGR entering the galaxy. 


Table 1: Magnetic held contributions to UHEGR scattering. The parameter 'i?p20 is for a 
proton with E = 10^° eV. The IGMF entry represents the smallest scattering angle possible 
with the parameters considered in this work. 
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2 kpc 
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Gal. Disk Tur. 

4 /iG 
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Gal. Halo Reg. 

2 ^G 

n/a 

8.7 kpc 

9.0° 

IGMF 

1 nG 

> 0.1 Mpc 

>4 Mpc 

>0.2° 


To give a sense of how the cosmic ray direction is affected by the galactic magnetic 
helds. Fig. 1 plots the mean scattering angle as function of the cosmic ray source energy for 
the halo regular, disk regular, and disk turbulent magnetic helds. The halo magnetic helds 
dominate the scattering while disk turbulent scattering tends to be a small effect. Iron nuclei 
are signihcantly dehected even for the galactic magnetic held disk turbulent contribution. 
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Galactic Magnetic Field Scattering 



Figure 1: Mean scattering angle vs. cosmic-ray energy plots for different galactic magnetic 
field contribntions (marked by line style). The thin lines are for protons while the bold 
lines are for iron nnclei. See Table 1 for the parameters defining each galactic magnetic field 
contribntion. See the text for a discnssion of the range of validity of these parameterizations 
and how they are nsed in this stndy. 


For the IGMF, we set Bq = 1 nG and provide resnlts for three coherence length valnes 
\b = 0.1, 1, 10 Mpc according to the constraints provided by [26] and [27]. Fig. 2 plots 
the mean scattering angle as fnnction of the energy for different source distances: 10 Mpc 
and 75 Mpc for protons and iron nuclei. Iron nuclei display significant deflection even for 
relatively nearby sources. 
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Figure 2: Mean scattering angle-energy plots for various combinations of magnetic field 
coherence length Xb, cosmic-ray charge Z, and propagation distance D. See the text for a 
discussion of the range of validity of these parameterizations and how they are used in this 
study. 


While the studies leading to the scattering models used here provide the average or root- 
mean-square scattering angles, they do not provide or discuss their statistical distribution. 
For regular helds, the scattering angles are deflected systematically by the mean scattering 
angle given by Eqn. 2.1. For turbulent helds, we assume that the scattering angles (Eqn. 
2.2) are Rayleigh distributed with mode (Tgcat, which is related to the mean scattering angle 
via 

( 2 . 

The scattering from a source catalogue is estimated by sampling this distribution. 



3 Cosmic Ray Source Model and Propagation 

We model the sources of UHECRs as emitters following a power law energy spectrum with 
constant cosmic ray luminosity. The arrival hux and energy at Earth are affected by inverse- 
squared distance losses and attenuation due to interactions with cosmic photon backgrounds. 
The following section describes the parametric model used to account for propagation energy 
losses. 
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3.1 Cosmic Ray Source Model 


We assume the cosmic ray source luminosity spectrum LgrciEsrc) (the rate of particles emitted 
at a given energy) for source energy Eg^c foi' each source in the catalogue follows a power law 
with index 7 ^ 

Lsrc{Esrc) = Lq {Esrc/E q) (3.1) 

where Lq is assumed to be a universal source luminosity at a reference energy of Eq = 10^® eV. 
We consider values of 7 ^ = 2.0 and 7 ^ = 2.7 but present results only for the latter. We found 
no significant difference in our results between the two values of 7 ^ considered. 

We assume a detector with effective area ^e//,det exposure time T that is the same 
for each source indexed by i. The total number of particles with an arrival energy at the 
observation point E^bs that is above a cut energy of E^ut is given by 


^ ^ /•oo 

E(^EQbs > ^ ^ / dEgrc^i^Eghgi^Egrc) E^^i) 

i=i 


Lsrc (E S 


(1 + Zi)47rd‘^{z., 


(3.2) 


In this equation M is the number of sources in the catalogue. The function EobsiEgrc) is 
the arrival energy resulting from a particle with energy Eg^c at the source at redshift Zj. 
The propagation is described in the subsections below. dc{zi) is the comoving distance of 
the source. The function Q{Eobs{Esrc) — Ecut) is the Heaviside step function requiring that 
the energy of arrival to the observer be greater than the cut value. See Appendix A for a 
derivation of Equation 3.2. 

For this study, we sample the catalogue to obtain the desired number of events assuming 
the detector has necessary exposure. The sampling is performed following the sum and integral 
in Equation 3.2. Each source in the catalogue is sampled according to [{1 + Zi)d‘^{zi)]~^, which 
in practice is achieved by rejection sampling. The source energy Eg^c is sampled according to 
the inverse power law with index 7 ^, and propagated to the energy of arrival to the observer 
Eobs- If Eobs > Ecut, the event is accepted, otherwise it is rejected. 


3.2 Cosmic Ray Energy Loss Length 

Cosmic ray energies are attenuated through particle interactions as they propagate through 
the Universe. Protons at the highest energies lose energy by producing pions through inter¬ 
actions with the CMB [1, 2]. The energy loss length for protons, shown here in Fig. 3, shows 
that UHECRs with energies above 6 x lO^^eV are likely to be substantially attenuated. We 
adopt the proton loss length calculations by [40], based on the analytical formulae from [41]. 

Heavy nuclei, such as iron, lose energy through pair production on the CMB and photo¬ 
disintegration through interactions with CMB and IR-UV background photons [42]. Photo- 
hadronic energy loss has been re-examined in recent years, as they are based on empirical 
measurements of the intergalactic background radiation and photonuclear interactions. We 
use [43] as our model for iron energy loss length, which is representative of the modern cal¬ 
culations [44-46]. 
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Figure 3: Energy loss length as fnnction of cosmic ray energy. This plot has been imple¬ 
mented by considering [40] for protons {Z = 1) and [43] for iron nnclei (Z = 26). 


3.3 Cosmic Ray Energy Propagation 

As the cosmic ray propagates throngh space, it is snbject to adiabatic losses dne to the 
expansion of the Universe as well as losses dne to interaction with backgronnd photons. The 
density and energy of backgronnd photons is also changing as the Universe expands. To 
acconnt for these losses, we discretely propagate the cosmic ray energy from its sonrce at 
redshift z to z = 0 in steps of Az. Given the energy Ej at step j, the energy Ej^i at step 
j -|- 1 is given by 

The second term in the right hand side is dne to the adiabatic redshift losses while the third 
term is given by the interaction with backgronnd photons corrected for their evolntion as 
the Universe expands. The energy-dependent cosmic ray attenuation length dne to photon 
interactions, discnssed in Section 3.2, is given by X^{E). The Hnbble constant today is given 
by Hq while the Hnbble parameter at redshift z is given by H{z). Equation 3.3 is derived in 
Appendix B. 

Fignres 4 and 5 plot the cosmic-ray energy and energy loss length as a fnnction of 
propagation distance assnming sonrces located at varions redshifts z. Assnming the energy 
loss length models for protons and iron nnclei in Fignre 3, it is evident that for sonrces with 
z > 0.03 cosmic rays reach Earth with energies below 6 x 10^® eV. Cosmic rays observed above 
this 6 X 10^® eV mnst originate from nearby sonrces. 

Fignres 4 and 5 demonstrate that the energy loss length thronghont the trajectory of a 
cosmic ray can vary signihcantly depending on its sonrce redshift. This occnrs becanse the 
backgronnd photon density and energy distribntion evolve with redshift (see Appendix B), 
which also means that the energy loss length (Fig. 3) evolves with redshift. For example, at 
sonrce redshift z = 0.3, the photon density and mean energy at the beginning of the cosmic 
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ray trajectory are lower than for a particle at source redshift z = 1.0, resulting in a longer 
energy loss length. As the cosmic rays propagate from sources at different redshifts, they 
will lose energy at different rates. The photon density and energy distribution along their 
trajectories also change at different rates for each case. 
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Figure 4: Examples of cosmic ray proton propagation for a source at redshift 2=0.01, 0.03, 
0.1, 0.3, and 1.0 at energies E = 10^^, 10^^, and 10^^ eV. The top panels show the proton 
energy evolution as a function of comoving distance of propagation. The bottom panel shows 
the energy loss length as a function of distance. This varies depending on the source redshift 
due to the changes in background photon density. 
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Figure 5: Examples of cosmic ray iron nuclei propagating for a source at redshift z=0.01, 
0.03, 0.1, 0.3, and 1.0 at energies E = 10^^, 10^^, and 10^° eV. 

In Fig. 6 we show a simulated distribution of cosmic ray energies at Earth obtained from 
sampling 10,000 events from an jg = —2.7 source luminosity spectrum. The distances have 
been sampled from the VCV catalogue assuming all sources have the same luminosity (see 
next Section for details on sampling). The distribution of energies for iron nuclei shows a 
dramatic cutoff for energies > 5 x 10^^ eV compared to protons. The location of the cutoff is 
due to the dominant contribution of nearby sources in the VCV catalogue, the nearest being 
at ~4 Mpc. 

4 Correlation to the VCV source catalogues 

Two approaches may be used in looking for correlations among known sources and UHECR 
arrival directions. One is to use a complete, uniform survey of matter in the universe. High 
correlations among such maps would indicate that UHECR production follows from regions 
of high star formation and gas density. The second approach is to look for correlations with 
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Energy Distribution of Protons 
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Figure 6: Energy distribution of cosmic rays propagated to Earth for 10,000 samples of 
sources sampled according to the VCV catalogue distance distribution assuming an 
luminosity energy spectrum. (Top) Proton UHECR distribution energy at the source and at 
the Earth. (Bottom) Iron UHECR distribution energy at the source and at the Earth. 


catalogues of candidate UHECR acceleration sites, such as active galactic nuclei (AGN) or 
gamma ray bursts (GRBs). 

The 2MASS survey [47] and the PSCz surveys [48] use infrared and near infrared ob¬ 
servations to make the most uniform maps available. Dust in older, dimmer galaxies and 
proto-galaxies alike emit in the infrared, so they are the best tracers of large scale structure. 
Similarly, Berlind and others [49] have generated simulation catalogues based almost entirely 
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on the matter density and ACDM N-body simulations. 

The Veron-Cetty & Veron (VCV) catalogue 12th edition [50] is a compendium of known 
AGN, largely derived from the 2dF catalogue and the Sloan Digital Sky Survey. While it is 
known to be a non-uniform survey, we use it here based on the precedent set by prior searches 
for correlations with AGN by Auger [3, 5, 6], HiRes [7], and TA [8]. However, hereafter, we 
treat it as a mock catalogue of UHEGR sources, by using it both as the source distribution 
and catalogue for correlation analysis. Our objective is not to test whether UHEGRs are 
produced by AGN in the VGV catalogue, but rather to characterize the ability to correlate 
given the effects of scattering and energy losses. 

4.1 Characterization of the VCV Catalogue 

Treating the VCV catalogue AGN as the sources of ultra-high energy cosmic rays, we charac¬ 
terize the expected angular separations between sources. The goal is to estimate the angular 
uncertainties required to possibly distinguish between catalogue sources. In this work, we will 
assume an optimistic angular resolution of 1°, given that orbiting fluorescence detectors such 
as JEM-EUSO [10, 51] expect angular resolutions between 1° and 3° and that radio-detection 
satellites [11, 12] expect angular resolutions of ~ 1°. We model the detector angular resolution 
by sampling a Rayleigh distribution around the arrival direction. 

The ability to distinguish between VCV catalogue AGN can be characterized by looking 
at the nearest neighbor angular distance of each AGN. In Fig. 7 we show the distribution of 
AGN nearest neighbor angular distance for both the catalogue with sources <75 Mpc and 
the inverse-square distance sampled version. Both distributions are similar and we find that 
~l/3 of sources have a neighbor within 2°. Within an angular distance of 5°, 70%-80% of 
sources have a neighbor in the VCV catalogue. This level of source confusion would make 
it difficult to identify individual AGN catalogue sources with high confidence, regardless of 
magnetic field scattering effects. 



AGN nearest neighbor angular distance [° ] 


Figure 7: Veron-Cetty & Veron catalogue 12th edition AGN nearest neighbor angular dis¬ 
tance before and after inverse-square distance sampling. 


Another important effect to take into account, assuming the source AGN in the catalogue 
have similar luminosities, is the inverse-square distance source rate reduction described in 
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Section 3.1. A distant source will contribute a lower flux on Earth compared to a nearby one 
simply due to the flux-luminosity relationship. It can be seen in Fig. 8 that the inverse-square 
distance sampling of VCV catalogue sources alone favors nearby sources. Despite the increased 
abundance of far away sources, the probability density that sources at a distance >100 Mpc 
contribute to the arrival flux on Earth is reduced by more than an order of magnitude and 
becomes negligible past 300 Mpc. 


Distribution of Distances 
for VCV Catalogue AGN 



AGN distance [Mpc] 


Figure 8 : Veron-Cetty & Veron catalogue 12th edition AGN distance distribution before 
and after inverse-square distance sampling of the catalogue. For this hgure, a source is only 
allowed to be sampled once. Without this requirement (as done in the simulations performed 
here), the contribution due to the nearest sources signihcantly outweighs the rest. 

It is worth mentioning that the number distribution of AGN as a function of distance in 
the VGV catalogue is not representative of the source evolution of AGN within a few hundred 
megaparsecs from Earth due to selection effects. The VGV sources, used as a source catalogue, 
would correspond to a strong negative redshift evolution that is unlike any ultra-high energy 
cosmic ray source candidate population. This approach leads to an increased contribution 
from nearby sources, with reduced intergalactic magnetic held scattering, compared to more 
distant sources. The lower bound resulting from this approach is therefore likely optimistic 
compared to using a source redshift evolution of a complete AGN catalogue. To avoid the 
complexities associated with densifying the VGV catalogue we use its listed sources as they 
are published keeping in mind that, although perhaps optimistic, this approach results in a 
perfectly valid lower bound. 

4.2 Source Correlation Analysis 

In this section, we describe the correlation of cosmic rays sampled from the VGV catalogue 
to VGV catalogue sources. The goal is to determine under which scenarios it is possible to 
discriminate between a catalogue correlation and an isotropic source distribution with > 5a 
conhdence as a function of cosmic ray sample size. 

The correlation analysis follows the procedure applied by the Auger collaboration [3, 5]. 
We compare the arrival direction of an event with the position of the source in the catalogue. 
An event correlates to the catalogue if the angular distance between the source and the 
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arrival direction of the cosmic ray is within a correlation angle, V'. We will explore the effect 
of varying this parameter in our simulations. Given a sample of N cosmic ray events, k of 
which correlate to the catalogue, the probability of the data being correlated to the catalogue 
is given by that fraction Pdata = k/N. We apply a cut on the distance of sources used for 
correlating against the event sample. We call this the correlation distance cutoff. Based on a 
data-driven optimization study, the Auger collaboration [3, 5] hxed the correlation angle at 
ijj = 3.1° and the correlation distance cutoff to 75 Mpc. The correlation distance cut restricts 
the sources used in the correlation analysis to 50% of the total number of sources in the VCV 
catalogue. 

The conhdence interval for Pdata is estimated from a binomial distribution P{pdata) = 
Cb{k, N)p^ata^^ —Pdata)^~^, where Cb{k, N) is the binomial coefficient. We calculate pdata for 
both a simulated arrival direction data set generated from the source catalogue and from an 
isotropic distribution. The isotropic distribution randomizes the arrival directions according 
to a uniform spherical distribution. We calculate the difference, Apdata, between pdata for the 
catalogue-generated distribution and for the isotropic distribution. The statistical signihcance 
from isotropy is then calculated as the number of conhdence intervals between pdata for the 
catalogue-generated distribution and for the isotropic distribution. We express the statistical 
signihcance in units of a, the Gaussian-distribution-equivalent number of standard deviations. 

The parameters tested in this study are source composition and the properties of the 
galactic and intergalactic magnetic helds. We look at the effects of galactic and intergalactic 
magnetic helds by testing the following scenarios: 

• GMF halo and disk (regular+turbulent) components and IGMF. 

• GMF disk (regular+turbulent) component and IGMF. 

• IGMF only. 

in addition to testing the IGMF Xb at 0.1, 1.0, and 10 Mpc for a magnetic held strength of 
Bq—1 nG in each case. We also compare source compositions comprising 100% protons and 
100 % iron nuclei. 

Rather than describing in detail the results of each scenario, we will hrst give an example 
of the analysis for source correlation analysis of 1000 proton events with arrival energies greater 
than 60 EeV detected with an angular resolution of 1 degree, assuming that the protons are 
scattered in both the GMF halo and disk. All of the scenarios considered here assume a 100% 
detector efficiency with respect to energy and perfect energy resolution. Our results place a 
lower bound on the sensitivity to source catalogue correlations. Results with other simulation 
parameter values of interest are shown in Tables 2 and 3. 

The probability of correlation, pdata, as a function of correlation angle, iIj, for both the 
VGV catalogue and an isotropic source distribution are shown in Fig. 9. It is clear that 
an optimal value of iIj exists that statistically discriminates between source catalogue and 
isotropic source distribution correlations. In Fig. 9 we have plotted Pdata vs. assuming 
different values of the IGMF coherence length Xb- We hnd that there is no significant effect 
on the optimal value of ip that depends on A^. 
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Figure 9: Probability of correlation as function of correlation angle ■0 for 1000-proton UHE- 
CRs, with 0.1, 1 and 10 Mpc IGMF coherence length, scattering in both the galactic halo 
and disk, and 1° detector angular resolution. The dashed line is for isotropic events and the 
solid line is for catalogue sampled events. The shaded regions are 5a confidence intervals. 


- 15 - 

















The optimal value of iIj is found by subtracting the Pdata for an isotropic source distri¬ 
bution from Pdata for the catalogue and finding the angle corresponding to the maximum, 
fpmax- Tables 2 and 3 report the optimal correlation angles for all scenarios considered in this 
study. The optimal correlation angle is below ten degrees in all cases, justifying our use of 
small scattering angle approximation in Eqns. 2.1 and 2.2. Regardless of the model used for 
large deflection angles, simulated events with ^ > 10 ° do not correlate to their source. 

4.3 Results 

Table 2: Statistical significance from isotropy, in units of a for selected proton and iron 
scenarios, for energy index 7 = 2.7, assuming a detector resolution of 1°. The last column 
contains the number of sigmas away from isotropic distribution for 100 , 1000 , 10000 proton 
UHECR events. 


z 

IGMF AsfMpc] 

GMF 

Resolution [°] 

'4^max [ ] 

N== 10 ^ 

N== 10 ^ 

N^^IO'^ 

1 

0.1 

Halo & Disk 

1 

6.36 

1.1 

4.0 

13.5 

1 

0.1 

Disk Only 

1 

4.14 

2.9 

11.2 

39.1 

1 

0.1 

None 

1 

2.78 

3.9 

17.4 

63.1 

1 

1.0 

Halo & Disk 

1 

6.41 

1.2 

4.1 

13.3 

1 

1.0 

Disk Only 

1 

4.26 

2.5 

10.4 

36.4 

1 

1.0 

None 

1 

3.35 

3.3 

13.9 

50.0 

1 

10.0 

Halo & Disk 

1 

6.50 

1.1 

3.9 

13.0 

1 

10.0 

Disk Only 

1 

5.03 

2.1 

7.7 

26.1 

1 

10.0 

None 

1 

4.79 

2.3 

8.2 

27.9 

26 

0.1 

Halo & Disk 

1 

8.30 

0.4 

1.0 

3.4 

26 

0.1 

Disk Only 

1 

7.90 

0.6 

2.0 

6.6 

26 

0.1 

None 

1 

5.92 

1.4 

5.4 

17.4 


Table 3: Statistical significance from isotropy, in units of a for selected proton and iron 
scenarios, for energy index 7 = 2.7, assuming a detector resolution of 3°. The last column 
contains the number of sigmas away from isotropic distribution for 100 , 1000 , 10000 proton 
UHECR events. 


Z 

IGMF AsfMpc] 

GMF 

Resolution [°] 

V’maxl ] 

N= 10 ^ 

N= 10 ^ 

N^^IO"^ 

1 

0.1 

Halo & Disk 

3 

6.51 

1.2 

4.0 

13.1 

1 

0.1 

Disk Only 

3 

4.95 

2.2 

8.4 

28.4 

1 

0.1 

None 

3 

4.64 

2.6 

9.2 

31.5 

1 

1.0 

Halo & Disk 

3 

6.53 

1.2 

3.8 

13.1 

1 

1.0 

Disk Only 

3 

5.04 

1.9 

8.0 

27.1 

1 

1.0 

None 

3 

4.75 

2.2 

00 

bo 

30.1 

1 

10.0 

Halo & Disk 

3 

6.58 

1.0 

3.9 

12.9 

1 

10.0 

Disk Only 

3 

5.48 

1.8 

6.7 

22.4 

1 

10.0 

None 

3 

5.36 

2.0 

6.8 

23.6 

26 

0.1 

Halo & Disk 

3 

8.30 

0.3 

0.9 

3.3 

26 

0.1 

Disk Only 

3 

7.90 

0.6 

2.0 

6.5 

26 

0.1 

None 

3 

6.07 

1.4 

4.8 

16.4 
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Tables 2 and 3 lists the significance of the correlation to the parent catalogne over an 
isotropic distribntion for several possible IGMF coherence lengths, galactic magnetic field 
models, and composition models. The significances are compnted assnming a 1° detector 
angnlar resolntion in Table 2 and a 3° detector angular resolution in Table 3. 

There are several proton-dominated scenarios which result in highly significant source 
catalogue correlations. Inclusion of galactic magnetic field scattering reduces the likelihood 
of a significant detection (> 5a) the most. On the order of 1,000 events are required for 
scenarios with only protons and full scattering off the galactic halo and disk to discriminate 
between anisotropic and isotropic distribution of UHECRs at the 5a level. In all proton 
cases considered, scattering in the halo substantially reduces the significance of a detection 
of correlations. 

There are no realistic scenarios wherein iron UHECRs result in significant correlations 
with the catalogue. We only report the results for Xb =0.1 Mpc, because the results for longer 
coherence lengths are insignificant except for 10,000 UHECRs that experience no scattering 
in the galactic disk or halo. Even the optimistic case that assumes no scattering within the 
galaxy, small (0.1 Mpc) ICMF coherence lengths, and a sample size of 1000 events shows 
correlations only at the 4 — Scr level. If scattering in the disk is included, then at least several 
thousand events are required to detect significant correlations. If scattering in the halo is 
also included, then no significant departure from isotropy is found for any number of events 
considered. 

The effects of varying the coherence length are most pronounced in scenarios where we do 
not include galactic magnetic field scattering. In such cases, it is clear that longer coherence 
lengths scatter both protons and iron more, making them more consistent with isotropy. For 
cases that include scattering in the galactic disk, the correlation significance of 1000 detected 
protons can be up to 4(T higher for Xb of 0.1 Mpc than for 10 Mpc. However, when including 
scattering in both the disk and the halo, galactic magnetic field scattering dominates over 
the ICMF contribution, such that the correlation significances vary by less than la for all 
coherence lengths. 

We find that an improved detector resolution does not translate directly to improved 
sensitivity to detecting a source catalogue correlation at greater than 5a. In the simplest 
cases without galactic magnetic field scattering, i/’max is reduced when the angular resolution 
is improved, thereby increasing the correlation significance. With scattering in the galactic 
disk, the improved angular resolution can reduce the required number of events to reach the 
5a level by a factor of a few. This effect degrades with increasing coherence length. However, 
when scattering in the halo is included, the differences in correlation significance are consistent 
with statistical fluctuations at the < la level. 


5 Conclusions 

The streamlined model presented here places a lower bound on the requirements for a state- 
of-art experiment to detect significant correlations from a source catalogue. We analyzed 
the required event rate above an energy threshold of 60 EeV of a future UHECR all-sky 
instrument for identifying sources, assuming several realistic scenarios with differing cosmic 
ray composition and magnetic field models. Such a simple parametric simulation does not 
require large computing power, but is capable of characterizing the trends and challenges for 
source identification. We have assumed perfect energy resolution, which neglects the effect of 
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cosmic rays below our energy threshold spilling over into the data sample, also making our 
lower bound optimistic. 

We hnd that when both the halo and disk magnetic helds are included in our scattering 
model, angular resolutions better than 3° do not greatly improve the detectability of the source 
catalogue correlations. The optimal correlation angle is greater than a few degrees in cases 
that include galactic magnetic held scattering. This is consistent with Eqns. 2.1 and 2.2 that 
scatter events with energies less than 100 EeV by several degrees. Combined with the power- 
law hux of UHECRs, most events will arrive at Earth scattered by more than 3°. Therefore, 
our results indicate that an improved angular resolution of 1° is not expected to signihcantly 
improve source catalogue correlations with our current understanding of scattering in the 
galactic halo. 

Scattering in the galactic disk isotropizes the cosmic-ray distribution more than scat¬ 
tering outside of the galaxy, despite the relatively unconstrained intergalactic magnetic held 
coherence length. Longer coherence lengths of the IGMF scatter UHECRs even further, but 
that scattering does not dominate the results of the correlation analysis. Future experiments 
would beneht from an improved understanding of the magnetic helds within the galaxy. This 
is underlined by how strongly the inclusion of scattering in the galactic halo affects the sig- 
nihcance of a correlation. 

If cosmic rays are predominantly iron, a detection of greater than lO'^ events above 
60 EeV would be required for source identihcation, and therefore, an exposure greater than 
100 times the state of the art. However, if they are predominately proton, an experiment that 
detects 10^ could expect correlations with a source catalogue at > 4cj even with dehection in 
the galactic halo and disk. This implies that a full sky survey of UHECRs should have an 
exposure at least 10 times the current state-of-the-art as well as an improved understanding 
of the composition of cosmic rays, which is consistent with the conclusions of [22]. Given the 
negative redshift distribution of the VCV catalogue, as used in this study, this is likely an 
optimistic lower bound. The simulations presented here indicate that the lack of a signihcant 
correlation of cosmic rays with energy > 60 EeV to nearby VCV catalogue sources in both the 
Auger and TA experiments does not exclude the possibility that AGN are the acceleration 
sites of UHECRs. 

A Number of Cosmic Rays Arriving from a Catalogue 

We calculate the number of cosmic rays observed at Earth from sources in a catalogue by 
assuming a specihc intensity at the source, total exposure on a given source, and tracking 
the energy losses due to propagation and cosmological expansion. We start from the specihc 
intensity /j, given as the number per unit time, area, solid angle, and energy, also known 
as differential hux, of cosmic rays arriving from source i to the observer. This is dehned 
according to the differential relation 

dN — li^fyifgdtobsdEobs di^fyifgdo'obsd^i^obs 

The dtobs term gives the arrival rate of UHECRs of energy between Fobs and Fobs + dFobs hr a 
detector differential area element daobs pointed at angle Bobs with respect to the source which 
subtends a differential solid angle di^i^obs in the sky. The specihc intensity is, in general, a 
function of the surface of the detector Fobs, and direction of observation Bobs, (pobs with respect 
to the source position. 
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The relation between the specific intensity of a sonrce as seen by the observer and that 
as seen from the sonrce frame of reference is given by 

sTc cos 0g'f'(^d(Tsvc d^^g 


^i,obs 


li 


(A.2) 


dtobs dEobs COS 0Q}^gd(7Qbs d^obs 

with the specific intensity of the sonrce illnminating the detector li^src defined by the relation 


dN = li^srcdtsrcdE, 


src COS ^i^src 


r.d^i 


(A.3) 


The nnmber of cosmic rays arriving at a detector with energy ifobs greater than Ecut 
from a catalogne with M sonrces is given by 


N{Eobs > Ecut) = 


M 

E 

i=l • 


dtobs I dEobs^i^^obs 


COS Oi^obsdCobsd^i ,obs di,obs 


(A.4) 

The sonrce time and observer time are related by dtsrc/dtobs = {^ + Zi), which gives 
the time dilation of the sonrce emission rate. The second term dEsrc/dEobs acconnts for 
the change in spectral band dne to energy propagation effects. The etendne of the sonrce 
dGsrc = cos OsrcdcTsrcd^src is the differential emission snrface area of the sonrce with a cosine 
projection factor on the solid angle snbtended by the detector, in the frame of reference of 
the sonrce. This is related to the etendne of the observer dGobs = cosOobsdctobsd^src via 
Etherington’s [52] reciprocity theorem dGobs = (1 + z/^dGsrc- We note that the distance 
between the sonrce and the observer from the point of view of the sonrce is the comoving 
distance, while from the point of view of the observer, it is the angnlar diameter distance. 

We may rewrite Eqnation A.4 as 


N{Eobs > Ecut) — 


M 


i + Zi 

1=1 


poo p p 

/ dEsrcQ{Eobs{E src ) Ecy^l ) / / 

Jo JAsrc J^i.i 


COS g/podcTg 


cdCli src Ii,i 


(A.5) 

where we have integrated over the observer time tobs to arrive at a total exposnre time T^. 
We have changed variables to Eobs, which is a fnnction of Egrc- The etendne is integrated 
over Asrc, which is the area of emission of the sonrce and over the solid angle snbtended by 
the detector as seen by the sonrce i^i^src- 

We assnme the specific intensity of each sonrce to follow the same model given by 


/ \ “T's 

Ii,src — df) f irsrc, dsrc, (/src) I 1 ■ (-^' 6 ) 

\ Efj J 

The scalar Iq sets the level of intensity. The fnnction f{rsrc, Gsrc, (/src) sets the emission beam 
pattern of the sonrce, which may, in general, depend on the location of the snrface of the 
sonrce Vgrc and direction of emission Ogre, (/src- We may write the specific intensity integrated 
over the sonrce etendne as 


cos Oi srcd(ysrcdCti^src di^sre — do 


E \ 

-d-'sre \ 


Er[ 


I j COS di^sredasred^i^sre f i^srci dgre-) 4^src^ 

J Asrc 

(A.7) 


to give 




cos Oj cyodcTcr'fd^^' 


2 ,src^^ sre^'^ ^t,src 


src — 


doA 


eff,src 


Aeff,det 

A7rdl,{Zi) 


Eo 


-"fg 


(A.8) 
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The solid angle subtended by the detector from the view of the source is i^i^src = 
Aeff^det/ ( 47 r(i^(zi)). The effective area of the source A^ff^src is the emission area referenced 
to an isotropic emitter (as is typically done in antenna theory). We write Lq = I^Aeff^src 
and assume the effective area of the source is the same for each source in the catalogue. We 
write the source luminosity as 


L'src{Egrc) — 


Eg. 


-la 


(A.9) 


to give 


M 


N{Eobs > Ecut) = 


Ti A 


eff,det 


^ 1 + Zj 47rd^(2;i) Jq 


dEgrc&iEgbgiE 

src ) ~ Ecut)L{E S 


(A.IO) 


Finally, we assume the exposure to each source is the same value given by TA(,ff^det and 
factor it out of the sum over all the catalogue sources to give Equation 3.2. 

B Energy Loss Propagator 

The energy loss propagator derivation follows the approach presented in [53]. The calculation 
has been adapted to make use of the energy loss length curves presented in Section 3.2. The 
energy losses associated with propagation has contributions from cosmological redshift as well 
as interaction with the photon background. Given that the photon background density also 
changes with redshift we model the energy change with redshift z as 


dE _ /dE\ 


+ 


dFy\ 


dz ydzj^^dshift ' \dzj^ 
where the energy losses due to redshift are given by 

(f-) ={l + z)-^E 

J redshift 

The losses due to photon background interactions ( 7 ) are given by 

dE\ _ /dE^ dx 
dzj \dx) dz 


(B.l) 


(B.2) 


where the energy loss is 




(B.3) 


(B.4) 


where is the loss length described in Section 3.2. The differential change in comoving 
distance dx/dz is given by 

“ (B.5) 


dz 


H{z)‘ 


The Hubble parameter is 


H{z) = Ho^J^Mil + zY + n^ 
for a ACDM cosmology with Hq—72 km/Mpc/s, VIm = 0.3, and Ha = 0.7. 


(B. 6 ) 
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The loss length also has cosmological corrections. The loss length X-yiE) for a cosmic 
ray with energy i? at z = 0 is given by 

1 /■“ da 

~ J {E,Ey) (B.7) 

Where -^^{E, E^) is the differential interaction cross-section of a cosmic ray of energy E with 
a photon of energy E^ and n^{E) is the photon density as 2 ; = 0. The loss length \'^{E) for 
a cosmic ray with energy E a. different redshift 2 ; is given by 


1 



dE'n'jE')%{E,E'). 


dE'^ 


(B. 8 ) 


The photon density at redshift 2 ; is given by n'{E!y) = {H{z)/Ho)^n{E^). Substitution of 
variables of E'^ for E^ in the integrand gives 



Combining the terms above into Equation B.l we arrive at 

dE 1 1 cHq 

dz 1 + z X^{E) H‘^{z) 


(B.9) 


(B.IO) 


We can discretize this equation for the energy Ej and redshift Zj at step j, the energy 
Ej+i at step j + ^ with redshift 2 ;j+i = Zj — A 2 ; is given by 


Ej+i — Ej 



A 2 ; cHq a 

x,{Ej)m{zj)) 


(B.ll) 
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